Tuning the electrical conductivity of nanotube-encapsulated metallocene wires 
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We analyze a new family of carbon nanotube-based molecular wires, formed by encapsulating metallocene 
molecules inside the nanotubes. Our simulations, that are based on a combination of non-equilibrium Green 
function techniques and density functional theory, indicate that these wires can be engineered to exhibit desirable 
magnetotransport effects for use in spintronics devices. The proposed structures should also be resilient to room- 
temperature fluctuations, and are expected to have a high yield. 
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The ability to fill carbon nanotubes (CNTs) with different 
compounds [1] has opened a new and growing field of re- 
search, where nanowire devices with new electronic proper- 
ties can be produced using simple chemical methods. Exam- 
ples include encapsulated nanowires made with metallic el- 
ements 01, with fullerenes that modify the local electronic 
structure of the nanotube 1 3], with metallofullerenes that mod- 
ulate the bandgap of semiconductor nanotubes and divide 
them into multiple quantum dots 0|, with organic molecules 
that allow to produce n- and p-doped nanotubes @], and 
with iron nanoparticles that facilitate the inclusion of mag- 
netic properties |7|,|8j]. These experiments demonstrate that 
it is feasible to fabricate long and stable chains of single 
molecules or atomic clusters by placing them in the interior 
of CNTs. They also show how the electronic and structural 
properties of the CNT can be conveniently tailored through 
the addition of different compounds. Recently, it has been 
demonstrated that metallocenes may be encapsulated inside 
single-walled CNTs with diameters of about 1 nanometer @|- 

Bis-cyclopentadienyl TM (TMCP2) is a metallocene com- 
posed of a transition metal atom TM, sandwiched between 
two aromatic rings made of five carbon and five hydrogen 
atoms each. It may be formed with the six elements in the 
middle of the 3d-row, V, Cr, Mn, Fe, Co, Ni, as well as with 
those in the column of iron, Ru and Os. All TMCp2 therefore 
share the same shape, which we depict in Fig. 1. Their elec- 
tronic structure is also very similar, with differences that can 
be traced to the diverse electronic filling of the d-shell of the 
TM atom. The crystal field associated with the benzene rings 
lifts the degeneracy of the d-shell into a triplet that lies lower 
in energy (d T 2_ y 2, d xy and d 3z 2_ r 2) and a doublet (d xz and 
dyz). The levels in the doublet are split into two peaks due to 
their bonding with the carbon atoms. In addition, there is an 
exchange field that breaks the spin degeneracy to accommo- 
date the magnetism predicted by Hund's first rule. Accord- 
ingly, VCp 2 , CrCp 2 , MnCp 2 , FeCp 2 , CoCp 2 and NiCp 2 have 
magnetic moments equal to 3, 2, 1,0, 1 and 2 Bohr magne- 
tons, respectively. FeCp 2 is paramagnetic since its core atom 
has a closed-shell electronic structure. The other metallocenes 
display a sort of mirror symmetry about FeCp 2 . To illustrate 
these points, we plot in Fig. 2 the density of states projected 
onto the d-shell of CoCp 2 . 



Such features are highly advantageous, since by placing 
metallocenes in the interior of CNTs, one may fabricate chains 
of molecules, whose electrical and magnetic properties can 
be tailored at will. We view the CNT as an external coating, 
whose conducting properties can also be modified by chang- 
ing its chirality from armchair to zigzag. Such a coating also 
provides a strongly enhanced chemical stability to the atomic 
chain and isolates it mechanically and, to a point, electrically. 
The simulations that we present here show that the formation 
energy of these wires is strongly exothermic for CNTs of radii 
between 4.5 and 7 A. We therefore argue that these nanomet- 
ric analogues of conventional electrical wires can readily be 
fabricated. A note of caution should be added however since 
fabrication and sorting of CNTs of a given chirality is cur- 
rently a very difficult task. 

In our search for metallocene nanowires, we have per- 
formed an extensive series of ab-initio simulations of TMCp 2 
inside CNTs of several chiralities. We have calculated their 
electronic and structural properties using the density func- 
tional theory code SIESTA 11 (ill , that uses pseudopotentials, 
and expands the wave functions of valence electrons by lin- 
ear combinations of numerical atomic orbitals. To be spe- 
cific, we have used a double-zeta basis set to describe hydro- 
gen and carbon atoms, and a double-zeta polarized basis set 
for the 3d metal atom. We have also used a number of k- 
points ranging between 20 and 100 to model the periodicity 
along the axis of the nanotube and the generalized gradient 
approximation to the exchange and correlation functional, as 
parametrized by Perdew, Burke and Ernzerhof fTll[l2Tl . We 
have relaxed the atomic coordinates of all atoms in the unit 
cell. Finally, we have calculated the electrical conductance of 
the devices with the non-equilibrium Green functions 
code SMEAGOL fl5ll . that has been successfully applied to 
the study of magnetorresistive effects in organic molecules 
and conductance quantization in atomic constrictions Jl6lfl7ll . 

As a first step, we have simulated single-wall CNTs of dif- 
ferent chiralities as well as isolated chains of TMCp 2 . We find 
lattice constants of 7.88, 7.99, 7.63, 7.55, 7.93 and 8.20 A for 
TM = V, Cr, Mn, Fe, Co and Ni, respectively. Metallocene 
chains can therefore be accommodated rather well inside arm- 
chair n) CNTs, whose theoretical lattice constant is 2.46 
A, with one molecule every three unit cells of the nanotube. 
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FIG. 1: (Color online) Schematic view of a TMCp 2 @ 3(7,7) nan- 
otube whith two encapsulated metallocenes, whose atoms have been 
highlighted in red, green and white colours (TM, C and H atoms, re- 
spectively). The axis of the farthest (closest) metallocene is aligned 
parallel (perpendicular) to the axis of the nanotube 

They can also be placed inside zigzag (n, 0) CNTs, that have 
a lattice constant of 4.32 A, with one molecule every two unit 
cells. Other chiralities may also be chosen for the CNT. For 
instance, (15,5) and (12,4) CNTs can also host metallocene 
chains (lattice constant of 15.59 A, in both cases), with two 
molecules per unit cell. 

To test the stability of the encapsulated chains, we have sim- 
ulated VCp2, FeCp2 and CoCp2 chains placed inside metallic 
(n, n) CNTs with n ranging from 6 to 10, which correspond to 
radii between 3.39 A and 6.78 A, respectively. The structural 
periodic unit of the simulations comprises N = 2, 3 or 4 unit 
cells of the nanotube and a single metallocene, and is denoted 
as TMCp2@^V(n, n). In addition, the axis of the metallocene 
has been oriented either parallel or perpendicular to the axis 
of the nanotube, as shown in Figs. 1 (a) and (b). We have 
calculated the binding energy Eb, per structural unit cell of 
the whole nanowire, defined as 

Eb = i?TMCp2@Af(n,n) — NE( n<n ) — i?TMCp 2 (1) 

as well as the charge transferred from the metallocene to the 
CNT AQ, and its magnetic moment m, as a function of the 
radii of the armchair CNTs. Fig. 3 shows the results for par- 
allel and perpendicular CoCp2@-/V(n, n) wires, with jV = 3 
and 4. We note that CoCp2@2(n, n) is energetically unfa- 
vorable since the metallocene molecules are too compressed, 
that CoCp2@3(n, n) accommodates rather well the metal- 
locene chain, and that the wave functions of the metallocene 
molecules have negligible overlaps in CoCp2 @4(n, n). 

Fig. 3 shows that parallel CoCp2@-/V(rt, n) nanowires are 
energetically favorable for radii equal or larger than about 4.5 
A. The most stable configuration is parallel CoCp2@iV(7,7), 
that corresponds to a radius of 4.75 A, and has a binding en- 
ergy as large as 0.85 eV. This result is in excellent agreement 
with the experimental data of Li and coworkers (j-J, who found 
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FIG. 2: (Color online) Densities of states gi of CoCp2, projected onto 
the d-shell orbitals of the cobalt atom for spin up (black solid line) 
and down (red dashed line). The index i denotes the orbital flavor, 

i — d x 2_y2 , dxz , C^3z 2 — r 2 ! dyz i dxy< 

a preferred radius of 4.67 A. The small energy differences be- 
tween N = i and N ~ A are due to the formation energy of 
the chain. 

We now examine the motion of single cobaltocenes along 
the nanotube axis. We note that large energy barriers will 
hinder the fabrication of chains with no faults, while if the 
barriers are too small, metallocenes will move almost freely. 
This will render the nanowires useless for operation at room 
temperature. We have therefore simulated several configura- 
tions of parallel CoCp2@4(n, n), where the cobaltocene is 
placed at different positions along the nanotube axis. We es- 
timate the energy barrier for motion of the cobaltocene Em 
as the energy -EcoCp 2 <»4(n,ri) of each configuration, referred 
to the energy of the most stable position. Thin nanowires 
have energy barriers of the order of several eV. On the con- 
trary, CoCp2@4(7,7) has an energy barrier of only 12 meV, 
which along with its large binding energy, makes this nan- 
otube a promising candidate for encapsulating cobaltocenes. 
For larger radii, Em decreases to values smaller than 3 meV, 
which allows the molecules to move almost freely along the 
nanotube, and eventually leave it OJ. 

To stablish wether or not metallocenes remain perfectly 
aligned at room temperature, we discuss now how easy is to 
rotate a cobaltocene from a parallel to a perpendicular con- 
figuration. We first note that perpendicular cobaltocenes only 
fit inside CNTs with n equal or larger than 8. We find that 
the energy cost to rotate the cobaltocene for CoCp2@3(8,8) 
nanowires is as large as 0.15 eV, and therefore they re- 
main aligned at room temperature. In contrast, metallocenes 
would be missaligned at room temperature in CoCp2@4(8,8) 
nanowires, due to the negligible overlap of metallocene wave 
functions that leads to a energy difference between both ori- 
entations of only 30 meV. 
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FIG. 3: (Color online) (a) Binding energy Eb, (b) charge trans- 
ferred to the nanotube AQ, and (c) magnetic moment m, of 
CoCp2 @N(n, n) wires, where n varies from 6 to 10. Blue 
solid and red dashed lines correspond to parallel and perpendicular 
CoCp2 @N(n, n), respectively. Circles and squares correspond to 
TV = 3 and 4, respectively. 



A Mulliken-population estimate of the total amount of 
charge transferred from the cobaltocene to the nanotube due 
to formation of the chemical bond, reveals that the nanotube 
is always n-doped as seen experimentally |9] and does not ac- 
quire any spin polarization. The amount of charge transferred 
is of about 0.3 to 0.4 e~, depending on the configuration, as 
shown in Fig. 3. The values of the magnetic moment are con- 
siderably smaller than those of the isolated molecule (1.0 /Xb), 
due to the charge transferred from the cobaltocene to the nan- 
otube @|. 

This charge transfer is also evident in the band structures 
shown in Fig. 4, where we plot the spin up and down elec- 
tronic bands of parallel TMCp 2 @3(7,7), with TM = V, Fe, 
Co. We note that the relative position of the CNT and the 
transition metal 3d bands indicate the amount of charge trans- 
ferred and therefore the strength of the chemical bond. The 
iron metallocene shows three 3d bands at the Fermi energy 
Ep, that are essentially unsplit and dispersionless whereas the 
CNT bands show a large dispersion. This implies that the cur- 
rent flowing through these nanowires must be carried only by 
the external coating. Something similar happens for the vana- 
dium wire, where the bands are now exchange-split, but still 
dispersionless and therefore should not conduct. In contrast, 
the 3d bands of the cobaltocene nanowire are exchange split 
and dispersive and should carry a spin-polarized current. 

The above discussion suggests that parallel CoCp2@3(7,7) 
and CoCp2@3(8,8) are the best candidates for spintronics de- 
vices and applications. On the one hand, these wires are 
the most stable devices from a structural point of view. On 
the other, their band structure shows the presence of a spin- 
polarized transport channel at the Fermi energy, which can 



give rise to a magnetoresistive effect. 

In order to test this possibility we have simulated devices 
comprised of two leads that sandwich a scattering region, and 
computed their conducting properties. The leads contain one 
unit cell of CoCp2@3(7,7) with spin up. The scattering region 
contains also one unit cell of CoCp2@3(7,7) with spin up or 
down. This makes a total of nine nanotube unit cells and 315 
atoms for the simulation. We denote each device by the spin 
(u or d) of the metallocene in the scattering region. Flipping 
such spin from up to down costs an energy of 0.05 eV, which 
implies that the ferromagnetic configuration can be stabilized 
at room temperature but that a relatively small magnetic field 
may switch it. We note that we have not included in our sim- 
ulations contacts between the devices and external electrodes, 
that should decrease to some extent their conducting and mag- 
netorresistive behavior. 

The u configuration, in which the spin of the three cobal- 
tocenes points upwards, corresponds to a perfect ferromag- 
netic chain, where there is no scattering. In this case the trans- 
mission coefficient is determined from the band structure and 
is equal to the number of open scattering channels. Fig. 5 
shows that there are three channels for spin up in a region 
around the Fermi energy that extends from -0. 1 eV up to 0.4 
eV. Two of them correspond to the nanotube and one to the 
cobaltocene. This excludes a tiny segment, where both cobal- 
tocene bands contribute. On the other hand, there are only two 
channels at the Fermi energy for spin down, that correspond 
to the nanotube, while the cobaltocene channels appear 0.3 eV 
above Ep. 

In order to study the change of resistance as a function of 
the magnetic configuration we have also calculated the con- 
ductance of the d configuration. In this case, the cobaltocene 
channel has moved to energies way above Ep due to the ex- 
change splitting, whereas the transmission for spin down is 
unaffected. As a consequence, the difference in conductance 
between u and the d configurations is 0.9 e 2 /h. Defining the 
magnetoresistance ratio as SG = (G u — Gd)/Gd, the pre- 



FIG. 5: (Color online) Conductance of the CoCp2@3(7,7) magne- 
torresistive devices for (a) spin up and (b) spin down electrons, u, 
d, dd and ddd conductances are plotted with black solid, red dashed, 
green dash-dotted and blue dash-double-dotted lines, respectively. 
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FIG. 6: Band structure of parallel CoCp 2 @2(11, 0), 

VCp 2 @2(ll,0) and 2CoCp 2 @l(i2, 4) for spin up (top pan- 
els) and down electrons (bottom panels). 



dieted ratio is SG = 20%. 

Since the spin flipping of only one cobaltocene may be dif- 
ficult to achieve in practice, we have simulated two other con- 
figurations that have two and three CoCp2 @3(7,7) unit cells in 
the scattering region, with all its metallocenes pointing down- 
wards. These configurations, that have a total of 420 and 525 
atoms, are denoted by dd and ddd, respectively. As shown in 
Fig. 5, the resulting transmission for spin up channels around 
the Fermi energy decreases slightly. This yields an additional 
increase in the magnetorresistance ratio to 25%. We have 
also simulated devices where the spin of the leads are either 
aligned or anti-aligned, but where there is no cobaltocene in 
the scattering region. In this case the magnetorresistance ratio 
vanishes. This confirms the fact that the spin polarized chan- 
nel is present only when the cobaltocenes form an unbroken 
conducting chain, with dispersive bands. 

Finally we have examined the effect of the conducting prop- 
erties of the CNT coating. With a view to improving the 
magnetoresistance ratio, it is desirable to increase the frac- 
tion of the current passing through the metallocene chain and 
so instead of using armchair CNTs, which are always metal- 
lic, the question arises of whether we could engineer an in- 
sulating coating by choosing non-conducting zigzag or chi- 
ral CNTs. Our calculations show that the charge transferred 
from the metallocene to the CNT due to the chemical bond- 
ing brings the CNT valence bands down in energy. The CNT 
coating becomes either conducting or semimetallic, as we 
show in Fig. 6 for CoCp 2 @2(ll, 0)), VCp 2 @2(ll,0) and 
2CoCp 2 @l(12,4). 

In summary, we have demonstrated that electrical 
nanowires consisting of a core made of a metallocene chain 
and a CNT coating may be fabricated, whose magnetic and 
electrical properties can be tailored by the appropriate choice 
of metallocene molecule and chirality of the CNT. We have 
performed a series of ab-initio simulations of the electronic 
and structural properties of several TMCp2@(n, n'), which 
agree very well with recent experiments and show that engi- 



neering of molecular spintronics properties is possible. 
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